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Abstract 

A brief review is presented of the stability domain of three- and four- 
charge ground-states when the constituent masses vary. Rigorous re- 
sults are presented, based on the scaling behavior and the convexity 
properties deduced from the variational principle. They are supple- 
mented by accurate numerical computations. 
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1 Three unit-charge systems 

Beyond hydrogenic atoms, the simplest systems in molecular physics consist 
of three unit charges, ±(+1, —1, — !)• The Hamiltonian reads 

if = ^ + ^ + ^-^-^ + ^. (1) 
2mi 2m2 2m3 ris 

A non-relativistic point-charge interaction is assumed, and everything else 
is neglected: relativistic corrections, finite-size effects, annihilation, strong 
interaction, etc. 

Nevertheless, as pointed out by Thirring it is not completely obvious 
why some systems are stable, e.g., = {p,p,e~), Ps^ = (e"'",e~,e~) or 
H~ = (p, e~,e~), while others, such as {p,p,e~) or (e~,p, e"^), immediately 
dissociate into a neutral atom and an isolated charge. 

One can use the scaling properties of (|1]) and similar Hamiltonians to 
choose ^ = = 1, as well as J^ca = where = l/rrii is the inverse 
mass of particle i. It is natural to use inverse masses, as the Hamiltonian 
depends linearly upon them. The condition cti = 1 is reminiscent of energy 
conservation in three-body decays, any possible configuration of which can be 
represented in a "Dalitz plot", which is an equilateral triangle of unit height. 
Each inverse mass is read as the distance to a side, as shown in Fig. [l|. 




Figure 1: The domain of normal- 
ized inverse masses a,. 




Figure 2: Schematic shape of the 
stability domain. 



Initially, some incorrect results where obtained for these systems, though 
these were sometimes based on very astute reasoning. It is, indeed, a more 
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delicate operation to estimate numerically where the stability frontier is than 
to calculate the binding energy of a system whose stability is well established. 
Recent numerical calculations are, however, very accurate and reliable. See, 
e.g, Refs. 0, 1^, H for some examples and the results in the literature, and 
Ref. P] for a comprehensive survey. If one summarizes the results of the 
literature, one finds that the stable configurations belong to a band around 
the symmetry axis m2 = where like-sign charges bear the same mass. 

The shape of the stability domain, as schematically pictured in Fig. |^, 
results from three basic properties 0: 

1. Any state along the symmetry axis is stable. This was shown by Hill 
0] using the trial wave function 

^ = exp (-ari2 - bris) + {a ^ 6} , (2) 

where the possibility of having a ^ b is crucial, as already noted 
by Chandrasekhar ^ for H~ and Hylleraas for Ps". Note that 
if a wave function with only an adjustable effective charge Ze, i.e., 
^ cx: exp(— Ze(ri2 + ^13)), or even more generally, any factorized wave 
function f{ri2)f{ri^) is used, it is not possible to demonstrate the sta- 
bility of H~ [|l^; one needs the anticorrelation a 7^ 6, or an explicit 
dependence, or both. 

2. Each instability region is convex. This results from the convexity prop- 
erties of the ground-state energy, after suitable rescaling such that the 
threshold energy Elmf ,1712) (or E{mi ,m^) on the other side) is kept 
constant. 

3. Each instability region has as a star shape with respect to A2 or A3. 
For instance, if one draws a line from A3 toward the symmetry axis, it 
corresponds, after rescaling, to increasing while the threshold energy 
E{mf ,7712) remains constant. Increasing decreases the three-body 
binding energy and thus improves stability. 

An application of this approach can be made when searching the critical 
masses for which stability disappears, for instance, 

a: (m+,oo",l") (3: (oo+,m^,l") 7: (l+,m^,l"). (3) 

In Fig. ^ one can see that a safe (lower) bound on rria and a safe (upper) 
bound on result in a straight line entirely in an instability domain. It 
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Ai 




As 

Figure 4: Schematic shape 
of the stabihty domain of 
{rrii, 1712 , ''^3 5 ''^4 ) states in the 
Figure 3: Some critical points tetrahedron of normahzed in- 

of interest. verse masses ai = m~^, with 

Etti = 1- 

intersects at 7' the mi = axis, this providing a upper bound for the 
critical mass m^, which turns out to be better than the limits obtained by 
direct studies of the particular configurations (l+,m~, 1~). In other words, 
global considerations on the triangle of stability are sometimes more powerful 
than local studies. 



2 Four unit-charge systems 

Consider now (mj*", m^, 7713 , 7714 ) systems, whose Hamiltonian reads 

+ ^ + i + + ^ + (4) 

2mi 2m2 2m^ 2m4 rig rss r^^ ' 

with, again, some scaling properties, so that c^i = 1 can be imposed: each 
configuration is a point in a regular tetrahedron, and the inverse constituent 
masses = represent the distances to the faces. 
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The lowest threshold consists of either two atoms or one three-body ion 
and an isolated charge. The lowest two-atom state is obtained by combining 
the heaviest of the positive charges mi and m2 with the heaviest of the 
negative charges and m^, and allowing the two lightest masses to form 
a second atom. This can be seen explicitly from the Bohr formula, but this 



is a more general results for any given potential, as shown in Refs. |Tl], |T2 . 
If for instance the first three particles are heavy and form a stable ion, then 
the lowest threshold is {mf ,m^) + mj. As the ion attracts the last 
charge with a potential that is asymptotically attractive and Coulombic, and 
thus supports many bound states, one is sure in such situation that the four- 
body system is stable. In practice, stability is guaranteed by finding a wave 
function whose expectation value for H is smaller than the energy of the 
lowest two-atom threshold. 

For the equal-mass case, i.e., the positronium molecule Ps2, this was done 
in 1947 by Hylleraas and Ore |]13[. Later on, it was realized that the stability 



of Ps2 implies that of all hydrogen-like configurations (M"*", M"*", m~, m~), 
VM/m. This includes the ordinary hydrogen molecule, starting from M = 
m, instead of the usual (and more natural) M — oo limit of the Born- 
Oppenheimer-Heitler-London treatment. The reasoning, made indepen- 



dently by Adamowski et al. |TJ] and in Ref. [0, is reminiscent of the well- 
known argument that if an even Hamiltonian, say /leven = p"^ + x"^, is sup- 
plemented by an odd term /lodd, then the ground state is lowered. This is 
seen explicitly if /lodd = x, but the general result follows simply from the 
variational principle applied to heven + ^odd using the even ground-state of 
heven as trial wave function. Here, parity is replaced by charge conjugation, 
and the (M"*", M"*", m~ , m~) Hamiltonian is rewritten as 



H — Hcvcn + Hodd — 
+ 



(5) 



The full Hamiltonian H has a lower ground-state than its even part Hcveni 
which is nothing but a rescaled version of the Ps2 case. The nice feature is 
that H and -ffeven have the same threshold, since when one computes the two- 
body energy, one first averages the inverse masses and to calculate 
the reduced mass. 

Once the binding energy of {M~^ , ,m~ ,m~) is given, then minimal 
extensions of the stability domain beyond the line {mi = m2, = m^} can 
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be derived |]T6|, ff? 



Many other configurations have been studied in the hterature. In particu- 
lar, all those with ms = m4 have been found to be stable. This property was 
checked extensively in Ref. |T^. Thus, for two identical electrons {m^ = m^), 
the molecule is always stable, irrespective of whether are the masses of the 
nuclei are equal or unequal, heavier or lighter than the electrons, or even one 
lighter and one heavier. 

The results obtained can be very conveniently described using the tetra- 
hedron of Fig. 



The thick lines 



which is also referred to by Armour 
correspond to stable configurations with ai = 02 = or 03 = ^4 = 0. The 
middle is just the familiar hydrogen molecule with infinitely heavy protons. 
The thin sides are outside the stability domain. The configuration at the 
middle of a thin line corresponds to {p,e~^ ,p,e~), studied in detail in |jl9|. 
The protonium {pp) forms a neutral and compact object with much too small 
an internuclear separation to bind the positronium atom (e"*"e~). 



3 Further studies 

So far, we have restricted ourselves to the ground state and to unit charge. 
Several generalizations can be envisaged. 

The case of three particles with arbitrary charges ±(gi, —q2, —qs) is dis- 
cussed rather extensively in Refs. |2^, |2l[| . 

In principle, our study can be repeated for excited states. For the first 
state with specific quantum numbers (negative parity, or a given value of the 
total angular momentum), the convexity properties remain, as one is dealing 
with the ground state in another sector of the Hilbert space. For the second 
or third state of given quantum numbers, the situation is more delicate. 

For excited states, stability should be addressed with respect to the ap- 
propriate threshold. For instance, a three-charge state of unnatural parity 
= 1+ cannot decay into the ground state of an atom and an isolated 
charge, unless a photon is also emitted. If one neglects radiative processes, 
then the lowest threshold involves an excited atom and is thus much higher 
than the threshold relevant for natural-parity states. 

The rigorous results on these systems are rather limited. There are for- 
tunately more and more powerful methods which can be used for numeri- 
cal investigations, among which are the stochastic variational method 



and Monte-Carlo methods [23]. Among the remarkable results obtained by 
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sophisticated numerical methods, we aheady mentioned the evidence that 
(m^,m^, 1~, 1~) systems are stable Vmi, m2- There is also the discovery of 
stable excited states for the Positronium molecule [p^ . 

Coulomb systems with more than four particles have been investigated. In 
the case of bosons, say (m"*", m~)^ , it has been known for long-time that the 
system becomes more and more bound when the number 2N of constituent 
increases. A kind of precocious approach to the large N limit was found in 
Ref. ||2^. More precise binding energies are provided in Ref. [^. In more 



recent papers, Varga et al. revealed several new stable configurations with 
an atom or an ion binding two positrons ||2^, ^ . 



The methods can also be applied to other potentials. For power-law 
potentials, the scaling laws are comparable to those of the Coulomb case, and 
the stability region can be drawn in a triangle for three-body systems, and 
a tetrahedron for four bodies, with similar convexity properties. Otherwise, 
the representation acquires additional dimensions. 



The role of the number of space dimensions has also been studied ||28 
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